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Summary. Shell effects in atomic nuclei are a quantum mechanical 
manifestation of the single-particle motion of the nucleons. They are directly 
related to the structure and fluctuations of the single-particle spectrum. Our 
understanding of these fluctuations and of their connections with the regular or 
chaotic nature of the nucleonic motion has greatly increased in the last decades. 
In the first part of these lectures these advances, based on random matrix theories 
and semiclassical methods, are briefly reviewed. Their consequences on the 
thermodynamic properties of Fermi gases and, in particular, on the masses of 
atomic nuclei are then presented. The structure and importance of shell effects in 
the nuclear masses with regular and chaotic nucleonic motion are analyzed 
theoretically, and the results are compared to experimental data. We clearly 
display experimental evidence of both types of motion, leboeuf@lptms.u-psud.fr 
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1 Introduction 



The mass is one of the most basic properties of an atomic nucleus. In spite 
of its fundamental character, it is a quite complicated and highly non-trivial 
quantity. According to Einstein's celebrated law, m = E/c 2 , all types of 
energy stored inside a nucleus contribute to its mass. If we imagine that 
the nucleons are initially separated, and then progressively move towards 
each other to finally form a nucleus in its ground state, a certain amount of 
energy B will be released in the process. Mathematically, this corresponds to 
K^ , the equation, 

U ■ M c 2 = Vm 3 c 2 -B , (1) 

1 

where M is the mass of the nucleus, c the speed of light, and rrij the mass 
of the jth nucleon. The binding energy B is responsible for the cohesion 
of the system, the larger B is, the more stable is the nucleus. There are 
different sources to this energy. The most important one comes from the 
strong attractive interaction between nucleons. Other interactions and effects 
that contribute to B are the Coulomb repulsion between protons, the surface 
effects, etc. In phenomenological descriptions of the mass, these contributions 
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are usually taken into account by, for example, liquid drop expressions a la von 
Weizsacker pQ. There is another important contribution. In a semiclassical 
picture where each nucleon is thought of as an individual particle moving in a 
mean-field potential, this contribution is related to the motion of the nucleons 
inside the nucleus. Because of Pauli's exclusion principle, when the nucleons 
are put together to form a bound system they are not at rest, but instead move 
with a speed which is of the order of 0.3c. Via Einstein's relation, this kinetic 
energy also contributes to B, and therefore to the mass. Part of this energy, 
namely the one that varies smoothly with the number of nucleons, is already 
taken into account in the phenomcnological liquid drop terms mentioned 
above. In contrast, the remaining part of the kinematic energy fluctuates with 
the number of nucleons. It moreover depends on the nature of the motion of 
the nucleons. For this reason, we refer to this contribution as the "dynamical" 
component of the mass. 

From the theory of dynamical systems, there are two important, extreme, 
and distinct possibilities for the motion of the nucleons inside the nucleus, 
namely it can be either regular or chaotic. A natural question to ask concerns 
the relation between these two types of motion and the properties of the 
associated dynamical mass. Is it possible to distinguish between a "regular" 
and a "chaotic" mass? How can the difference be detected experimentally, if 
any? The purpose of these lectures is to bring some answers in this direction. 
We will show that it is indeed possible to make the difference between both 
types of dynamical masses via their fluctuation properties, i.e., how the mass 
varies as a parameter, typically the number of nucleons, is changed. When the 
results are applied to the available experimental data, the analysis of both 
types of dynamical masses indicates that, apart from a dominant regular 
component that will be described in some detail via a spheroidal potential, 
there might be in addition some chaos in the ground state of the nucleus. The 
corresponding contribution to the mass and its dependence with the number 
of nucleons will be explicitly computed. It will be shown that the hypothesis 
of chaotic layers explains the observed differences between the experimental 
data and previous theoretical calculations. 

Although the setting we are using is based on a single-particle picture of 
the nucleus, several arguments indicate that the results obtained are of more 
general validity. They may in fact be related to correlations acting between 
nucleons that are beyond a mean-field picture. 

Traditionally, dynamical effects in the structure of nuclei are referred to as 
shell effects. This designation originated in atomic physics, where the symme- 
tries of the Hamiltonian produce strong degeneracies of the electronic levels 
(the shells). These degeneracies induce, in turn, oscillations in the electronic 
binding energy. Shell effects are therefore due to deviations (or bunching) of 
the single-particle levels with respect to their average properties. The degen- 
eracies of the electronic levels of an atom produced by the rotational symme- 
try are an extreme manifestation of level bunching. In general, in systems that 
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have other symmetries, or have no symmetries at all, there will still be level 
bunching, but its importance will typically be minor. Therefore, depending 
on the presence or absence of symmetries, the shell effects may be more or 
less important. The level bunching, and more generally the fluctuations of the 
single particle energy levels, are thus a very general phenomenon. The theo- 
ries that describe those fluctuations make a neat distinction between systems 
with different underlying classical dynamics (e.g., regular or chaotic). It is 
therefore our purpose, before entering in the discussion of the nuclear masses 
in 21 to give a general (though elementary) presentation of the fluctuation 
properties of the level density (|2|and|sJJ, and of the way these fluctuations 
manifest in the different thermodynamic properties of a Fermi gas (2J- It will 
be shown, in particular, that generically the size of the fluctuations is more 
important when the classical underlying motion is regular, as compared to 
the chaotic case. This fact may seem paradoxical, since we are usually willing 
to associate chaos with noise and fluctuations. However, concerning the shell 
effects in quantum mechanics, the situation is exactly the other way around: 
because of the instability of the classical orbits, level clustering is less impor- 
tant in chaotic systems. We will precisely quantify the difference for several 
thermodynamic quantities. We will moreover discuss universality, e.g. the 
validity in large classes of systems of common features of the fluctuations. 

The "dynamical" or shell fluctuations of the mass [21 El considered here 
are a particular example of a general phenomenon. Similar fluctuations are 
expected to occur, with different degrees of importance, in all thermodynamic 
quantities of a fermionic gas. This point is discussed in some detail in 21 
Many illustrative examples may be mentioned, like the fluctuations of the 
persistent currents in mesoscopic rings, the force fluctuations observed when 
pulling a metallic nanowire, or the supershell structure in metallic particles. 
They all have the same physical origin, which in the context of semiclassical 
theories is associated to classical periodic orbits. See for example Ref. |3] for 
a discussion in condensed matter mesoscopic physics, and Ref. (5] for cluster 
physics. For a general statistical theory of the fluctuations, and further related 
references, see Ref. 0. 

These lectures are based on research carried out over the past few years, 
during which I have enjoyed and benefitted from collaborations with A. 
Monastra and O. Bohigas. I am deeply indebted to them. 21 is based on 
Leboeuf and Monastra, Ref. |JJ, and 21 0n Bohigas and Leboeuf, Ref. [7J. The 
analysis of the regular component of the mass and of supershell structures 
in nuclei presented in 21 based on a spheroidal cavity with a finite lifetime 
for quasiparticles, is original and was not published elsewhere. References jHJ 
and [HI report some closely related results that are not discussed here; they 
treat thermodynamic aspects of a fictitious element, "The Riemannium", a 
schematic many body system inspired from number theory. 
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2 Local Fluctuations: Random Sequences 

Consider a bound single-particle Hamiltonian H , whose quantum mechanical 
spectrum is given by a discrete sequence of energy levels {Ei, E%, . . . , Ei, . . .}. 
H may either represent a self-consistent mean field approximation of an 
interacting system, or simply the Hamiltonian of a given one body problem. 
At a given energy E, we denote the typical distance between neighboring 
eigenvalues, the mean level spacing, by 5. The aim in this section is to present 
a short overview of some of the results that have been obtained concerning 
the fluctuations of the single-particle energy levels on scales of order 8, and 
to establish connections with the underlying classical dynamics. Fluctuations 
on a scale S are termed "local" , compared to fluctuations on larger scales to 
be discussed in J3 

One of the most pervasive theories in the description of the statistical 
properties of sequences of numbers is the random matrix theory (RMT). Its 
range of applicability largely exceeds the spectra of single-particle systems, 
and even the frontiers of physics. For example, one can find it in the descrip- 
tion of nuclear, atomic and molecular systems, the motion of electrons in a 
disordered potential, the behavior of classically chaotic systems, the study 
of integrable models, the description of the statistical properties of the crit- 
ical zeros of the Riemann zeta function, and in problems of combinatorics. 
Some interesting articles and lectures covering these topics can be found in 

Refs. nm-di]. 

Motivated by its mathematical simplicity, one of the original ensembles in- 
troduced is the so-called Gaussian ensemble of random matrices [Rj\ , defined 
as the set of hermitian matrices H whose elements are Gaussian indepen- 
dent variables. The probability density of a given N x N hermitian matrix is 
defined as 

P(H) = c H &cp(-^trH 2 Y (2) 

where cr is a normalization constant. A very specific aspect of this ensemble 
is its invariance under rotations in Hilbert space. The form of H depends on 
symmetry considerations: its matrix elements are real for even spin systems 
with time reversal symmetry, complex in the absence of time reversal sym- 
metry, and quaternion for odd spins with time reversal symmetry. For these 
three cases, the parameter fi in Eq. (J2J takes the value 1, 2 and 4, respectively. 
Given the probability density of matrix elements J2J), the problem is to 
find the probability distribution of the associated eigenvalues Ei, i — 1, . . . , N 
and eigenvectors {pk} of H (the number of parameters in the parametrization 
Pk of the eigenvectors depends on the symmetry of H). In the basis where 
H is diagonal, trH 2 = ^E 2 . Moreover, the Jacobian of the transformation 

i 

from matrix elements to eigenvalues and eigenvectors is given by |16| 
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The eigenvector's components are absent in this expression because of the 
rotational invariance of the ensembles. From these results, it follows that the 
joint probability density of the eigenvalues is (in an abuse of notation, we use 
the same symbol P as for the matrix probability density) , 
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(4) 



where ce is a normalization constant. 

The most characteristic feature of this probability density is the existence 
of a strong repulsion between eigenvalues. The repulsion may be seen, in 
Eq. J2J, from the fact that P — > when Ej — > 2?*. Being related to the 
Jacobian, this repulsion is not a particular feature of the ensemble l|2*jl. but 
rather constitutes a very basic and general fact of most quantum mechanical 
systems. 

From Eq. J3J the fc-point correlation function, 



R k {E 1 ,...,E k )=N\/(N-k)l 



dE k +i ■ ■ ■ dEjyP(Ei, . . . , En) , 



may be explicitly computed. We are particularly interested in the case k = 
2, that defines the density of pairs of eigenvalues separated by a distance 
y = Ei — E\ (when the distance y is local (i.e., on a scale 5), the two- 
point function is stationary and depends only on the difference; moreover, we 
are here interested in the "bulk" results, i.e. statistics of eigenvalues located 
near the center of the spectrum, far from the edges). A convenient way to 
characterize the function Ri(y) is through its Fourier transform, 



K{t) = 2ft / cos(yr/ft) [R 2 (y) -f+p 6(y)] dy , 



(5) 



usually called the form factor. K{t) has units of time; p is the average density 
of states, 

P = 5- 1 =T H /h, (6) 

where we have introduced the Heisenberg time, the conjugate time to the 
mean level spacing 6. From Eq. J3J, the random matrix form factor is found 
to be P| . 



Krratir) = < 



2r - r log (l + £)] 0(t h -t) + 
2r H -rlog(^A]9(T-T H ) 



k T0(t h -t)+t h 9{t-t h ) 



[3=1 



(3 = 2 



(7) 
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For simplicity, we have restricted to (3 = 1 and 2 (since the spin dynam- 
ics will not be included in the following, we do not consider the symplectic 
symmetry). In 0) the function is Heaviside's step function. t h is the only 
characteristic time scale in Eq. 01. For short times r -C t h , the form factor 

behaves as, 

2 

K rmt {T) = - T , T « T H , (8) 

whereas for r 3> r H it tends to t h (this is a general property of any discrete 
sequence of levels) . 

Another useful quantity is the nearest-neighbor spacing distribution p(s), 
defined as the probability to find two neighboring eigenvalues separated by 
a distance s, where s is here measured in units of the mean level spacing 6. 
Though the limiting behavior N — > oo of p(s) may be computed analytically, 
the result is not explicit and it is customary to use the expression obtained 
from 2x2 matrices (which turns out to differ, by accident, by only a few 
percent from the exact expression), 

ff S exp(-f s 2 ) 0=1 

p(s) = { (9) 

lf S 2 exp(-^) /? = 2. 

The results obtained from RMT are usually compared to an individual 
system, whose spectrum has been computed numerically or measured ex- 
perimentally. The statistical properties of the system are calculated in some 
window of size AE located around an energy E. To obtain statistically sig- 
nificant results, an average is done by changing the position of the window. 
The window size AE should be sufficiently large in order to include many 
single-particle levels and make a statistical analysis meaningful, but also it 
should not be too large, to guarantee that the gross features of the spectrum 
(e.g., the mean level spacing) are constant. 

The RMT, as formulated above, is a one parameter theory. The results 
(see Eq. 0) depend only on the mean level spacing S (or Heisenberg time t h ). 
In order to compare different spectra coming from different physical systems, 
the mean spacing should be normalized. Doing this, good agreement with 
RMT of the local fluctuations has been found in many different situations. 
This is the case for example of the neutron resonances in heavy nuclei, which 
was the original experimental motivation to develop the RMT. Much simpler 
though non-trivial problems were subsequently analyzed. At the beginning of 
the eighties it was shown that, in the metallic regime, the statistical properties 
of the eigenvalues of the famous Anderson problem (e.g., an electron moving 
in a metal with impurities, modeled by a random 3D potential) coincide with 
RMT 17 . Up to now, this is the only case were the agreement with RMT was 
proved rigorously, though the model is not fully deterministic (randomness is 
incorporated by hand in the Hamiltonian) . In fully deterministic systems, RM 
statistics were conjectured to hold generically for high-lying eigenvalues of 



Regularity and chaos in the nuclear masses 7 

systems with a classical fully chaotic motion. This is the Bohigas-Giannoni- 
Schmit (BGS) conjecture in quantum chaos ^3 \I§\, supported by many 
numerical results and by some analytic arguments (see next section) . One of 
the simplest Hamiltonians allowing to test the conjecture, and that will be 
here employed as a model system, is to consider the motion of a free particle 
inside a cavity with perfect elastic reflections on the boundary. The upper 
part of Fig.^shows several fully chaotic two-dimensional (2D) examples: the 
Sinai billiard, the stadium, and the cardioid. The quantization of the motion 
is given, for a billiard system, by the Schrodinger equation of a free particle 
with Dirichlct boundary conditions imposed (i.e., the wave function should 
vanish on the border). In the left part of Fig.|2|is shown the nearest neighbor 
spacing distribution of high-lying eigenvalues of the cardioid billiard obtained 
numerically The cardioid billiard is time-reversal symmetric, and is therefore 
compared with the ensemble (3 = 1 in Eq. © |2U| . 



Chaotic 
Motion 



Sinai 



Stadium 



Cardioid 




Circle 



Rectangle 




Ellipse 




Regular 
Motion 



Fig. 1. Examples of classical motion inside 2D cavities. Upper part: chaotic, lower 
part: integrable. In each cavity, a typical trajectory is shown. Taken from |21|. 



A relevant question is the universality of the ensemble of matrices con- 
sidered. There is no special reason, aside technical ones, to consider the par- 
ticular Gaussian distribution of matrix elements described above. Other sta- 
tistical weights could have been chosen as well. It is therefore important to 
understand to which extent the results are independent of the distribution. 
Some form of universality of the local statistics has been shown. As the size 
of the matrix increases, the bulk eigenvalue distribution converges towards a 
limit that, to a large extent, is independent of the probability measure used 

Through the BGS conjecture, RMT is associated to the concept of com- 
plexity, here manifested as chaotic classical motion. In contrast, there are 
large classes of physical systems were the conjecture does not apply. An ex- 
treme and well known class is that of the regular, or integrable, systems. 
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Fig. 2. Nearest neighbor spacing distribution for the cardioid and circle billiards, 
compared to RMT (/? = 1) and Poisson distributions. Taken from [21]. 



These are systems with as many independent integrals of motion as degrees 
of freedom. The lower part of Fig. ^ shows several two-dimensional (2D) ex- 
amples. In the circular billiard, for instance, the two independent integrals 
of motion are the energy and the angular momentum. The existence of these 
constants of the motion forces the orbits to form regular patterns, and no 
exponential sensitivity to perturbations, like in chaotic systems, exists. Reg- 
ular and fully chaotic motion are two extreme situations. A billiard with an 
arbitrary shape will generically contain both types of dynamics, regular and 
chaotic (mixed system). Some aspects of mixed dynamics will be discussed 
in |2 

When a regular system is quantized, the statistical analysis of the fluc- 
tuations give results that are quite different with respect to RMT. What is 
found is that the local statistics of the corresponding quantum energy levels 
are described by a Poisson process, (e.g., the eigenvalues behave, statistically, 
as a set of uncorrelated levels) |23 • Many numerical simulations support this 
conjecture, though there is no general proof, aside some particular systems 
25 . This property is illustrated on the right hand side of Fig. |2 where the 
nearest neighbor spacing distribution of a circular 2D billiard is compared to 
the Poisson result, p(s) = e~ s (6 = 1). The form factor of an uncorrelated 
sequence of levels is independent of r, 



l^unc\J~ ) 



(10) 



In summary, and loosely speaking because there are exceptions that we 
don't discuss here, for the two types of dynamics considered the local corre- 
lations of energy levels are universal and either described by random matrix 
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theory in chaotic or "complex" systems, or by an uncorrelated Poisson se- 
quence for integrable or "regular" systems. 



3 Long Range Fluctuations: Semiclassics 

In fact, and contrary to what Fig. [21 suggests, RMT or a Poisson process 
do not describe faithfully the fluctuation properties of real physical systems. 
The reason is that these simple ensembles fail to describe some important 
physical features related to the short time dynamics. There is only one time 
scale that can be associated, through quantum mechanics, to a stationary 
sequence of random numbers. This is the Heisenberg time, that corresponds 
to the conjugate time to the mean level spacing. But physical system possess 
other time scales, that are typically smaller than t h . In ballistic systems (e.g., 
no disorder present), like a particle moving freely inside a 2D or 3D cavity, 
another important time scale is the time of flight across the system, denoted 
Train- We will give a precise definition of r m i n below using the periodic orbits. 
This is typically the smallest relevant time scale of a dynamical system. It 
corresponds to a large energy scale, given by, 

E c = h/r min . (II) 

We will now see that the typical single-particle energy levels of a ballistic 
dynamical system have, superimposed to the universal local fluctuations on 
a scale 5 discussed in the previous section, system-specific long-range mod- 
ulations on a scale E c . These modulations of the density with respect to the 
average properties of the spectrum are not described by random sequences 
of numbers (correlated or not). 

It turns out that a very convenient tool to analyze in a systematic way the 
deviations of a quantum spectrum with respect to its average properties is 
the semiclassical analysis. This method is convenient from either a technical 
as well as a conceptual point of view, and is justified if a classical limit of the 
quantum problem exists. This limit can be obtained by different procedures. 
For example, in the context of nuclear physics it is known that a classical 
limit of some simple shell models of the nucleus can be obtained by letting 
the number of particles tend to infinity 26 27 . More generally, the connexion 
with a classical dynamics can be done through the mean field, which reduces 
the fully interacting many-body problem to a single-particle motion in a 
self-consistent potential. The interest of semiclassical theories is that they 
offer a unified scheme of analysis, and provide a conceptual framework that 
allows in many cases to understand a physical property using simple ideas. 

The basic object we are interested in is the single-particle density of states 



p(E,x)=g s J2S[E-E J (x)} . (12) 
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The Ej(x) are the discrete eigenvalues of the single-particle Hamiltonian H. 
The eigenvalues depend on a set of parameters x that fix, for example, the 
shape of the potential, or may represent any other external parameter. The 
prefactor g s accounts for spin degeneracy. 

The classical motion is described by a function H (p, q) of the phase- 
space variables (p,q). In the semiclassical limit fi — > the quantum density 
of states can be approximated by a sum of smooth plus oscillatory terms, 

p(E,x)=g s [p(E,x)+p(E,x)} . (13) 

The first term is a smooth function of E that describes the average properties 
of the spectrum. It is usually known as the Thomas-Fermi contribution (or 
Weyl expansion for the motion of a particle inside a hard-wall cavity) . The 
leading term of p is given by the well-known semiclassical rule that associates, 
in a _D-dimensional system, a phase-space volume hP to each quantum state, 

P=J^Djd D P d D q5(E-H(p,q)). (14) 

Corrections to Eq. (|14|) depend on derivatives of H [28] • In the case of the 
motion inside a 31? cavity, the expansion takes a particularly simple form 

129, 

p -^{-W) E 16^+16^ VW) E/ +°i E ')-^) 

Here V and S are the volume and the surface of the cavity, and £ is a 
typical length that depends on its topology; m is the mass of the particle. 
The integral of the density with respect to the energy is also of interest, 

Af(E,x)=g s J2&[E-E j (x)] ; (16) 

3 

it gives the number of single-particle levels with energy Ej < E (counting 
function). It can be decomposed as in Eq. 113|) . 

Af(E,x)=g s p7(E,x)+N'(E,x)\ . (17) 

The average part, obtained for a 3D cavity by integrating l|15|) with respect 
to E, is 

-r-, V /2m\ 3/2 , /9 S 2m C /2m\ 1/2 , /9 „, , 

7J=7r-s\-r) E/ -T7. t e + -e-\-t) E^ 2 + 0(1). 18 

6tt 2 V fi 2 J 16tt h 2 8tt V h 2 J \ J \ J 

Deviations with respect to the smooth behavior of the density are de- 
scribed by the fluctuating part p in Eq. I|13(l . They are given, to leading 
order in an ^-expansion, by |3()l I31j , 
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oo 

p(E,x) = 2^2^2A p , r (E, x)cos[rS p (E, x)/h + v Ptr ] . (19) 

p r— 1 

The sum is over all the periodic orbits p of the classical Hamiltonian H (p, q). 
The index r takes into account multiple traversals (repetitions) of a primitive 
orbit p. The orbits are characterized by their action S p — § p.dq, period 
t p = dSp/dE, stability amplitude A p ^ r , and Maslov index u pr . The functional 
form of A p _ r depends on the nature of the dynamics. In intcgrable system, 
the periodic orbits form continuous families, with all members of a family 
having the same properties. In contrast, in fully chaotic systems all orbits 
are unstable and isolated. This important difference is at the origin of the 
fact, to be discussed in the next sections, that shell corrections are more 
important in integrable systems than in chaotic ones. 

JNotei Why a sum over periodic orbits? This can be understood from the following simplified 
argument. The density of states may be calculated from the Green function, p(E) — — ^-Im I d q 
rim e ^ -|_ G(q, q, E-\-ie). In the Fcynman representation, G(q, q, E) is the propagator for paths of 
energy E that start at q and come back to q. In the limit h — > 0, the leading contribution to G is 
a sum over all the classical trajectories that start and come back to q. Finally, a stationary phase 
approximation of the above integral over q selects, among all the classical closed trajectories, 
those that start and come back to a given point with the same momentum. Those arc the periodic 
orbits. 

The fluctuating part TV of the counting function may be obtained by 
integration of Eq. I|19fl with respect to the energy. To leading order in a 
semiclassical expansion, the integration with respect to the rapidly varying 
factor Sp(E,x)/h dominates. This gives, to leading order, 

oo . 

Af(E,x)=2n^2Y2-^-sm[rSp/n + u p , r ] ■ (20) 

P r=i rTp 

To illustrate the semiclassical approximation of the density of states and 
of the counting function in a concrete example, we consider the free mo- 
tion of a spin-1/2 particle inside a three-dimensional (3D) cavity, elastically 
reflected off the surface. The spherical cavity and the corresponding semiclas- 
sical approximation have been extensively studied in the past, and I refer the 
reader to some of those references for a more detailed presentation |301 125) . 
Quantum mechanically, the single-particle energy levels E n j of a free par- 
ticle of mass m moving inside a cavity of radius R with Dirichlet boundary 
conditions are given by the quantization condition, 

= , n= 1,2,... , I = 0,1,2,... , (21) 

where E n ^ = h 2 k^ ;/2m = h 2 e„ t i/2mR 2 = E e n j, J is the Bessel function, 
and k the wavenumber. The previous expressions define the quantities e n j — 
(k n jR) 2 and Eq — h 2 /2mR 2 . The index I labels the angular momentum of 
the particle, and n, that serves to identified the different zeros of Jz+1/2 (£n,z) 
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Fig. 3. Counting function of a spin-1/2 particle in a 3D spherical cavity. Full stair- 
case: quantum mechanical. Dotted staircase: semiclassical (almost indistinguishable 
from the quantum mechanical one). Smooth interpolating curve: Weyl's formula. 
Detail in the inset. 



at fixed I, is the radial principal quantum number. The degeneracy of each 
quantum mechanical energy level is g n ^ = 21 + 1. 
The counting function is defined as, 



J^sph(e) = g s V" fjn.i 0(e- e n ,i) 



(22) 



which increases by a factor g s xg n ,i at each eigenvalue e ni ;. The exact staircase 
function is represented in Fig.[3](gs = 2), with the energy measured in units 
e = E JEq. Big jumps of the function at energy levels with large degeneracies 
are clearly visible. 

We now turn to a semiclassical description of Af sp h(e)- A classification 
and discussion of the periodic orbits of the sphere may be found in |30|. 
Each orbit is characterized by two integers, w (the winding number), and v 
(the number of vertices or bounces of an orbit off the surface) (see Fig. 0J . 
The semiclassical theory leads to an approximate expression of the counting 
function in terms of these orbits, given by Eqs. (|17|l . I|18|) and (|20|l . For the 
sphere, the smooth and oscillatory parts read |3(J| . 



A/" S ph(e) 



,3/2 
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,1/2 



9tt 



3tt 



-T-,. 



N sp h(e) = -2\/e }_^ 2_^ -j^- cos (yfe £ vw 

w=l v=2w 

The different factors entering Eq. (|24(1 are defined as, 

-(2™)- 1 



,tt/2) 



v = 2w 



(-l) to sin(2 vw )y/sm(6 vw )/irv e 1 / 4 u > 2w 



(23) 



(24) 



(25) 
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w + 3/2 
L,„, = 2v cos(6>„„ 



v = 2w 



v > 2w 



(26) 



I (27) 

V w = ttw/v . (28) 

The parameter £ vw is the length of the periodic orbit (v, w) of a unit radius 
sphere, and v vw is the phase attached to each orbit. The semiclassical ap- 
proximation of Af S ph is also represented in Fig. [3] as a dotted line, together 
with the average part Af sp h- Viewed on large scales, the accuracy of Weyl's 
formula is quite good. The inset shows a closer view of the function, where 
the fluctuations with respect to the average and the accuracy of the semi- 
classical approximation can be appreciated. Orbits with w up to 40 and v up 
to 80 have been used in Eq. IJ2"4l to obtain the curve in the figure. It is seen 
that the semiclassical approximation of J\f sp h is quite good. 






(2,1) 



(3,1) 



(4,1) 






(5,2) 



(7,2) 



(7.3) 



Fig. 4. Some of the periodic orbits of the sphere, labeled by (v, w), the number of 
vertices and the number of turns around the center, respectively. 



We now let for the moment the spherical cavity and return to general 
considerations concerning the trace formula. Equation i|19H shows that the 
discrete quantum mechanical spectrum is, formally, recovered from an in- 
terferent sum over all the periodic orbits whose contributions add up con- 
structively at the positions of the eigenvalues Ej, and destructively elsewhere. 
Since in the vicinity of some reference energy E* the action may be expressed 
as S P (E) m S P (E*) + t p (E — E*), each periodic orbit contributes to the den- 
sity with an oscillatory term, as a function of E, whose wavelength is h/r p . 
Long periodic orbits are therefore responsible for the small scale structures 
of the spectrum. To resolve the spectrum on an energy scale of the order 
5, and describe the departures from the average properties on that scale, 
long periodic orbits of period t h are needed. In contrast, short orbits de- 
scribe oscillations of the density of states on large scales. The shortest orbit, 
of period T m ; n , defines the outer energy scale E c = h/r m i n of the density 
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fluctuations. The density modulations on a scale E c are usually referred to 
as shell effects. Their importance heavily depends on the dynamical proper- 
ties and symmetries of the system. For instance, in systems with a central 
potential, the angular momentum conservation produce strong degeneracies 
among eigenvalues. The degeneracies create energy regions where the density 
of eigenvalues is high compared to the mean, and others where it is small, 
thus producing long range fluctuations of the density. This is an extreme il- 
lustration of level bunching. From a classical point of view, the enhancement 
produced by symmetries is due to the appearance of families of periodic orbits 
all having, because of the symmetry, the same action, amplitude, etc. This 
produces, at the level of the trace formula, an enhancement of the shell ef- 
fects. However, shell effects are present in arbitrary systems, with or without 
symmetries, since they are a manifestation of the existence of short periodic 
orbits, a generic property of any dynamical system. In less symmetric systems 
they will be less important because the degeneracy of the families of peri- 
odic orbits will be lower, or inexistent, like in chaotic systems where periodic 
orbits are isolated. Correspondingly, the associated shell effects will be less 
dramatic. This point will be reconsidered in the next sections. See Ref. [jjj- 

The scale E c is typically much larger than 5. A precise characterization 
of the difference is given by the adimensional parameter, 

g = E c /S = T H /r min , (29) 

which counts the number of single-particle energy levels contained in a shell. 
It is a measure of the collectivity of the long range spectral modulations 
imposed by the shortest periodic orbits. It will play an important role when 
quantifying the difference of size of shell effects in regular and chaotic systems. 
Denoting L — V 1 ^ 3 the typical size of the system, it follows from Weyl's 
law for a D-dimensional cavity that g is proportional to (kL) D ~ l , where k 
is the wavenumber at energy E. Thus, in the semiclassical regime of short 
wavelength compared to the system size, fei>l, the parameter g is large. 

To have a more precise idea of the consequences of the presence of an 
outer scale in the oscillatory structure of the density, and of its relation with 
the universal properties discussed in ^2 consider again the form factor. Using 
the semiclassical approximation of the density, Eq. Ijl9|l , it can be shown |32| 
that K{t) is expressed as, 



K{t) =h 2 (Y,J2 A ^ A P ', r > cos f ^ ''' ,S '"' 



\p-,p' r,r' 



{rr p +t j t p >) 



1 E 

(30) 
where the brackets indicate an energy average (see M4.2I for a more precise 
characterization of the averaging window) . 

For times r <C t h the off-diagonal contributions p ^ p' in Eq. I|30(l are 
eliminated by the averaging procedure, and the behavior of K(r) is well 
described by the diagonal terms p — p'. This gives a series of delta peaks at 
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' p; 



K{t) 



^E 

p.r 



//■-'> '< r *(T 



(31) 



The lowest peak is located at r = T mi „, and for r < r mi „ the form factor is 
identically zero. Since the position and amplitude of this peaks are system 
dependent, it follows that for times of the order of T min no statistical universal 
behavior of the form factor exists, and Eqs. © and (|10|l do not provide a good 
description. In particular, the fact that K{r) — for r < r m i n is completely 
out of reach of the models based on random sequences. 

Going to longer times, but still remaining in the regime of validity of the 
diagonal approximation, r m i n <C r <C t h , it can be shown that indeed the 
RMT behavior © in chaotic systems, and the uncorrelated statistics <|1(J[1 
in integrable systems, are recovered. Their semiclassical origin is related to 
the statistical behavior of long classical periodic orbits 33, 32 . There is no 
proof of the agreement for longer times, aside some additional corrections in 
the t expansion of the form factor in time-reversal chaotic systems |jil], and 
results in some particular integrable models (e.g., the rectangular billiard, see 
|25)). We will therefore simply assume that for times r ^> r min the statistical 
description of the previous section, based on random sequences, holds. Figure 
El summarizes the main features of the form factor for integrable and chaotic 
systems. 



K(T) 



Regular 



Chaotic 




Fig. 5. Schematic view of the form factor for regular and chaotic systems 



In summary: the fluctuations of the single-particle spectrum in ballistic 
systems are schematically of two types. On scales of order S, they are universal 
(described by an uncorrelated sequence for integrable systems, and by RMT 
in chaotic ones), as discussed in the previous section. In contrast, on scales of 
order E c = gS, with j> 1, there are long range modulations whose structure 
and amplitude are system specific. 
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4 Fermi Gas 

We now turn to the central theme of these lectures, the theory of the quantum 
fluctuations in interacting Fermionic systems. Our reference model will be a 
confined Fermi gas, i.e. a gas of non-interacting particles moving in a self- 
consistent mean field. It turns out, as it will be discussed below, that such 
an approximation is well justified when studying the thermodynamic fluctu- 
ations of interacting systems like the atomic nucleus or the valence electrons 
in metals. 

The starting point is a grand-canonical description, appropriate to study 
the properties of an open system. Although this choice is in part motivated 
by technical reasons, we will specify how to adapt the results to other phys- 
ical situations. In particular, to systems where the number of particles is 
conserved, like the atomic nucleus. 

When the chemical potential p, is fixed, the thermodynamic behavior of 
a Fermi gas is described by the grand potential, 

Q(fj,, T, x) =~k B Tf dE p{E, x) log[l + e (p-E)/k B T^ 

(32) 
= -JdEAf(E,x)f(E,n,T), 

where p(E, x) and J\f(E, x) are the single-particle density of states and the 
counting function, respectively, already introduced in JH T is the tempera- 
ture, and k B is Boltzmann's constant, x parametrizes the mean field potential, 
or represents some additional external parameter, and / is the Fermi function 

f^^ T )= 1 + J-, )/kB T - (33) 

When the number of particles A in the gas is fixed, the energy of the 

system is 

U(A, T,x) = JdEE p(E, x) f(E, fx, T) , (34) 

where p is determined from the equation 

A = j 'dE p(E,x) f(E,n,T). (35) 

Knowing these two thermodynamic potentials, other thermodynamic 
functions are directly obtained by differentiation with respect to p, x, or 
T. 

4.1 Semiclassical approximation 

We now exploit the semiclassical expansion of the density of states to obtain 
semiclassical expressions for the thermodynamic functions of the gas. Using 
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Eq. (|13|) , the density of states in Eq. 1(32(1 is written as a sum of smooth plus 
oscillatory terms. This approximation is valid in the semiclassical regime A ^> 
1 where the typical wavelength of the fermions at energy /i is much smaller 
than the system size, and is less accurate at the bottom of the spectrum. 

The substitution in Eqs. 1(32(1 and (|34|l of the density of states by its 
average behavior, Eq. 1)14(1 or l|15Jl . gives rise to the smooth, bulk classical 
expressions for the thermodynamics of the Fermi gas 35 . In finite fermionic 
systems, like the atomic nucleus or the electron gas in a metallic particle, it 
is however well known that the smooth part of the thermodynamic proper- 
ties of the interacting system are not correctly described by the Fermi gas 
approximation. In contrast, and this is a crucial point, it is known that the 
fluctuations of the properties of the interacting system with respect to the 
average behavior and, in particular, the energy fluctuations, are well approx- 
imated by the fluctuations of the Fermi gas. This is the content of the Struti- 
nsky energy theorem |3SJ. This property justifies, in interacting many-body 
Fermionic systems, the validity of the analysis of the fluctuations presented 
here, based on a single-particle mean-field picture. 

Equation 1)19(1 is then inserted in Eq. 1(32(1 to compute the oscillatory 
part of the grand potential. To leading order in h and for low temperatures 
{k B T <C p, degenerate gas approximation) the integral gives |S], 

fib, T, x) = 2h 2 g s J2 E r T^ cos ( r ^ + ^) ■ ( 36 ) 

p r—l P 

In this expression, all the classical orbits have energy E = p, and they all de- 
pend on the external parameter x. The temperature introduces the prefactor 

k t (t), 

M T ) = ■ ih x > r T = h/(2ir 2 k B T) . (37) 

Sinn (t/t t ) 

This factor acts as an exponential cut off for orbits with period t p > t t . For 
temperatures such that r m j„ 3> t t , the quantum fluctuations are washed out 
and only the smooth part of £2 given by p survives. 

Similarly we may analyze the energy of the gas at a fixed number of 
particles. The main technical difficulty here is that p is not constant as some 
external parameter is varied (like the shape of the nucleus for example) , but 
varies in order to satisfy, at any x, the condition 1(35(1 (p is in fact a function 
of A, T, and x). It should therefore also be decomposed into smooth and 
oscillatory parts, 

H = -p + Ji , (38) 

with 71 defined by the condition 

A=fdEp(E,x)f(E,JI,T). (39) 

The fluctuating part of the energy of the system is defined as 
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U = JdEE p{E) /(£, M , T) - JdEE p(E) f(E, % T) . (40) 

Using Eq. I|13[) . assuming that Ji <C ~p, and integrating by parts Eq. 1)40(1. one 
finds, 

U = JdEE p(E) f^(E,-p,T) jt-fdE Sf(E)J(E,JI,T) 

-JdEEAf(E)f E (E,JL,T), [ ' 

where ^ (E,Ji,T) = df(E,Jl,T)/d]X and f E {E,JL,T) = df(E,JI,T)/dE. From 
Eqs. 1|35(1 and (|39|) . in a similar approximation, it follows that 

/" d£ £ p(£) j' fl {E,Tl,T)p = JdEE M{E) f' E (E, JZ, T) . 

This implies, 

Z7(A, T, x) = - / e£E A/"(-E, x) f(E, Jl, T) . (42) 

Comparing this equation with Eq. (I32|) . we conclude that the fluctuations of 
the energy at a fixed number of particles A coincide, to leading order in an 
expansion in terms of Jl, with the fluctuations of fl evaluated at a chemical 
potential Ji(A 7 T, x) defined by inverting Eq. i|3"5|) . From Eq. l|3l)|) . we therefore 
get, 

U(A, T, x) = 2h 2 g s J2 E ^ 7 {l 2 ^ ™ S (rS p /h + v p , r ) . (43) 



T T 



All the classical orbits that enter this expression are at energy p(A, T, x). 

U are the variations of the energy due to kinematic effects (" dynamical 
fluctuations" ) . Equation (|43)l is the main theoretical tool to analyse its prop- 
erties. It associates to each periodic orbit a certain amount of energy that 
depends on its period, stability, etc. The total contribution is the sum of all 
the energies associated to the periodic orbits. The temperature enters only 
through the prefactor n T , whose magnitude decreases as T increases. For 
studying the contribution of kinematic effects to the mass of the nuclei, we 
simply take T = 0. Then k t — 1, and Eq. (|43|l becomes, 






A 



B(A, x) = 2h'g a ^ ^ ^ cos(rS p /H + ^ r ) . (44) 

p r—l P 

This equation describes the fluctuations of the mass due to dynamical ef- 
fects; they are added on top of the smooth part given by the von Weizsacker 
formula. At T = 0, fi is the Fermi energy E F of the system, which is also 
decomposed into a smooth and an oscillatory part, E F = E F + E F . The def- 
inition of the smooth part of the Fermi energy is, from Eq. 139JI at T = 0, 
given by 
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/E p pE ' £■ 

dE p(E, x)= dE p(E, x) = W(E F , x) , (45) 

which defines, by inversion, the function E F (A,x). All the classical orbits in 
Eq. (|4*4"j) are evaluated at energy E F (A, x) . 

The oscillatory part of other thermodynamic functions may be computed 
by direct differentiation of fl or U with respect to the appropriate parameter 

u 

4.2 Statistical Analysis of the Fluctuations 

Equation (|44|l requires the knowledge of the periodic orbits and of their dif- 
ferent properties (period, stability, Maslov index, etc) in order to compute 
the fluctuating part of the mass. In this Section, instead of a detailed com- 
putation and description of the fluctuations for a particular system, based on 
a precise knowledge of some set of periodic orbits, the aim is to study their 
statistical properties. The interest of such an approach is well known: using 
a minimum amount of information, a statistical analysis allows to establish 
a classification scheme among the fluctuations of different physical systems. 
It also allows to distinguish the generic from the specific, and provides a 
powerful predictive tool in complex systems. A general description of the sta- 
tistical properties of the quantum fluctuations of thermodynamic functions 
of integrable and chaotic ballistic Fermi gases, and of their temperature de- 
pendence, was developed in Ref. [Sj. Here we will only concentrate on the 
energy of the gas at zerojemperature, e.g., its mass. 

The fluctuating part B(A, x) shows, as a function of the number of parti- 
cles or, alternatively, the Fermi energy E F , oscillations described by Eq. (|44|l . 
The statistical properties of B, and in particular its probability distribution, 
are computed in a given interval of size AE F around E F . This interval must 
satisfy two conditions. It must be sufficiently small in order that all the 
classical properties of the system remain almost constant. This is fulfilled if 
AE F <C E F . Moreover, it must contain a sufficiently large number of oscilla- 
tions to guarantee the convergence of the statistics. As stated previously the 
largest scale associated to the oscillations is E c . Then clearly we must have 
AE F ^ E c . In the semiclassical regime the hierarchical ordering between the 
different scales is therefore 

8 < E c < AE F < E F . 

Second Moment and Universality 

The average value of the fluctuating part B defined as in Eq. I|44|) is zero. This 
is not strictly true in general, because it is possible to show that subleading 
corrections in the expansion l|42(l have a non-zero average. Ignoring this, 
the variance is the more basic aspect of the probability distribution of the 
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fluctuations. It provides the typical size of the oscillations, and can easily be 
compared with experiments. We will now compute a general expression for 
the second moment of the probability distribution, that allows also to make 
an analysis of the universal properties of the fluctuations. The key point is to 
understand which orbits give the dominant contribution in Eq. I|44|l . Though 
it is clear that the weight of an individual orbit decreases with its period, the 
net result of the sum of the contributions is unclear, because the number of 
periodic orbits of a given period grows with the period. 

From Eq. I|44|l the square of B is expressed as a double sum over the 
periodic orbits involving the product of two cosine. The latter product may 
be expressed as one half the sum of the cosine of the sum and that of the 
difference of the actions. The average over the term containing the sum of 
the actions vanishes, due to its rapid oscillations on a scale AE F . Therefore, 
letting aside for the moment the spin factor (which will be discussed later 
on), 

(fi 2 )=2^(^^^co S (^)\ . (46) 

\p,p' r,r' P P V 'I e f .x 

To simplify the notation we have included the Maslov indices in the definition 
of the action. Ordering the orbits by their period, and taking into account the 
restrictions imposed by the averaging procedure, the variance Eq. (|46|l can 
be related to the semiclassical definition of the form factor K(t), Eq. (|3L)[) . 
The variance of the mass fluctuations takes the simple form [7] , 



h 2 f°° dr 



<P') = ^l ? *M- ,47) 

Analogous expressions connecting the variance of different thermodynamic 
quantities, like for example the response of the gas to a perturbation, the 
entropy, etc, at any temperature, can be found in Ref. |H]. 

To obtain Eq. l|47fl we made use of the fact that the orbits giving a non- 
zero contribution to (|46[) have similar actions (unless their average will be 
zero). This implies that their period is also similar, and can be considered 
to be the same in the prefactor (but not in the argument of the oscillating 
part) . 

Based in Eq. I|47|) we now make a simple analysis of the variance of the 
mass fluctuations. From Eqs. 10 and IJBJ it follows that for chaotic systems 
the integrand in Eq. 1|47(1 behaves as t~ 3 for short times and t~ 4 for long 
times, while Eq. (|10|l implies that for integrable motion the integrand varies 
as r~ 4 . Therefore in both cases the integral (|47|l converges for long times. 
The dominant contributions come from short times, where the integrand is 
large. If the form factor of a pure random sequence is used in Eq. I|47|l . the 
integral diverges. In real systems the divergence of the integral is in fact 
stopped by the cutoff at r = T m i n of the form factor. Because, as shown in 
JU the short-time structure of the form factor is specific to each system (i.e., 
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it is not universal) , we see from Eq. (|47|l that in general the second moment 
of the mass fluctuations is non-universal, and consequently the same is true 
for the probability distribution. 

We thus conclude that the fluctuations of the mass are, regardless of the 
regular or chaotic nature of the motion, dominated by the short non-universal 
periodic orbits of the system. 

It should be mentioned that this is a particular property of the mass, 
that is non generic, and therefore not shared by all other thermodynamic 
functions. It can be shown that the fluctuations of some functions, like the 
entropy for example, are dominated at low temperatures by times much larger 
than T m i n ; the corresponding distributions are universal. For a given function, 
the universality can also depend on the nature of the dynamics. See Ref. |2]. 

4.3 The second moment in the T m ;„ approximation 

We have demonstrated that the main contributions to the dynamical part of 
the mass come from the shortest periodic orbits. On the other hand, we have 
seen in fJJJthat for short times a good approximation to the form factor (|30|l 
is to keep only the diagonal terms in the double sum over the periodic orbits, 
Eq. 1)31(1 . Using the latter approximation, the variance of the fluctuations of 
the mass is expressed as, 

(^-^E^. (48) 



P.!' 



V 



This expression requires the explicit knowledge of the periodic orbits. This 
information is not always available, like for example in the case of the atomic 
nucleus. There is, however, a simple way to estimate the variance that re- 
quires a minimum amount of information on the orbits. The approximation 
consists in using in Eq. (|47|l the corresponding form factor (of an uncorrelated 
sequence given by Eq. (|10f> for a regular dynamics, of RMT given by Eq. (7J 
for a chaotic motion), and to impose the additional and important condition 
K{j) = for r < T m i n in all cases. This is clearly an approximation, that we 
call the r m i„-approximation. It extrapolates the statistical behavior of the 
orbits down to times r ~ T m i n , ignoring the short-time system-dependent 
structures. All the short-time structures are condensed into a single param- 
eter, the period of the shortest orbit. 

The virtue of the r m i„-approximation is to provide a simple estimate of 
the size of the fluctuations, as well as of its dependence with the number of 
particles, using a minimum amount of information. 

Since the integral obtained from Eq. (|47|l in the r m i„-approximation is 
straightforward, we do not give here a detailed account of its computation. 
The result, for integrable and chaotic systems, given as an expansion in terms 
of the small parameter 1/g, is, 
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E 2 C (I- J +0(g- 2 log g)) Chaotic = 1 (49) 

RF^ (l - 3^) Chaotic = 2. 



It is evident from these expressions that the fundamental energy scale that 
determines the energy fluctuations is E C1 and not (5, or E F . This is natural, 
since we have shown that the energy or mass fluctuations are controlled by 
the long range fluctuations of the single-particle spectrum produced by the 
short orbits on a scale E c . Fluctuations of the levels in smaller energy scales, 
whose statistical properties are universal as discussed in SJ3J contribute with 
\/g corrections in Eq. H49JI . In chaotic systems the variance of the fluctuations 
is twice smaller in systems without time reversal symmetry. Semiclassically, 
this is also easy to understand, because in systems with time reversal symme- 
try each orbit is doubly degenerate (the primitive one and its time-reversed) , 
in contrast to systems with no time reversal symmetry. The coherent con- 
tribution of these pairs of orbits produces a variance twice larger for j3 = 1. 
Finally, the variance of the fluctuations is g-times larger in integrable systems 
compared to chaotic ones. This amplification, which could be quite large (a 
precise estimate for nuclei is given in the next section), has as semiclassical 
origin the existence of families of periodic orbits all contributing in phase in 
Eq. (|43|l . This estimate is valid for generic integrable systems. Even larger 
amplifications could exist in "super-integrable" systems (cf next section) . 

Higher moments of the probability distribution of the energy fluctuations 
may be computed similarly. Starting from Eq. 14411 . a generalization of the 
diagonal approximation used for the second moment can be implemented |H]. 
The results show that the moments of the distribution are generically all 
different from zero, giving rise to asymmetric probability distributions. For 
example, the third moment is given by, 

oo oc 
p n = lr 2 = l 

(50) 
where A p . r = A p , r /r 2 T 2 . This and the corresponding expressions for higher 
moments were explicitly tested in some integrable and chaotic models, were 
it is found that they work extremely well [S| |S1 03 . 



5 Nuclear Masses 

The binding energy of a nucleus is defined by Eq. flfl . It is a direct measure 
of the cohesion and stability of a nucleus. It is customary, as was done in the 
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previous sections, to analyze the experimental values observed by decompos- 
ing B into two parts, 

B{E,x)=B{E,x)-B(E,x) . (51) 

This splitting is at the basis of the so-called shell correction method, intro- 
duced in nuclear physics by Strutinsky [2Zj (see also 3, 38] ). Note the minus 
sign we have introduced in front of B. The average part of B is described by 
a liquid drop model a la von Weizsacker, 

— Z 2 (N — Z) 2 t 

B = a v A-a s A 2/3 - a c ^^ - a A ^ —^- - a p -^ . (52) 

The different terms in this expansion are associated to volume effects, sur- 
face effects, Coulomb interaction, proton-neutron asymmetries, and pairing 
effects, respectively. N and Z are the number of neutrons and protons, 
A = N + Z is the total nucleon number, and t\ — —1,0,-1-1 for even- 
even, odd-even, and odd-odd nuclei, respectively. Note that all the terms 
are smooth functions of A, N, and Z, except t\ that takes into account odd- 
even effects. 

In contrast to B, the oscillating energy B, that describes deviations with 
respect to the smooth part, can be analyzed, to leading order in a one- 
body density expansion, in terms of the fluctuations computed from a single- 
particle spectrum. This is the content of the Strutinsky energy theorem |36| . 
As explained in 21 t ms contribution is related to kinematic aspects of the 
nucleus. It contains explicit information about the motion of the nucleons 
in the self-consistent mean field. The net output of the semiclassical theory 
for the oscillatory part of the binding energy is Eq. (|44|> . that expresses B 
as a sum of contributions, each depending on a classical periodic orbit of 
the mean-field potential. Our purpose now is to explicitly investigate the 
ability of such a formula to describe the experimentally measured masses of 
the atomic nuclei. 

In order to proceed, the first thing to do is to subtract from the exper- 
imental data the average part B. If this is done from the 1888 nuclei with 
Z,N > 8 of the 1995 Audi-Wapstra compilation [HHI , the result is Fig. El 
The difference B = B — B is plotted as a function of the neutron number (top 
part). Each point in the figure represents a nucleus. The parametrization of B 
used is taken from Ref. |30|: a., = 15.67, a s = 17.23, a A = 23.29, a c = 0.714, 
and a p = 11.2 (all in MeV). 

A clear oscillatory structure (the shells) of amplitude ~ 15 MeV is ob- 
served, with well defined sharp minima, that contrast with the more smooth 
and rounded shape of the maxima. Another peculiar feature of the shells 
is that the wavelength of the oscillations clearly increases with N. A min- 
imum in this plot corresponds to a minimum of B (and therefore, locally, 
to a maximum of B). It defines more stable nuclei. The minima are located 
at the "magic" numbers N = 28,50,82, and 126 (N = 8 and 20 are not 
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Fig. 6. The mass difference B = B - B (in MeV) for the 1888 nuclei of the 1995 
Audi-Wapstra compilation. Top: as a function of the neutron number N. Bottom: 
as a function of iV 1 ' 3 . 



clearly visible in this figure). The usual interpretation of these minima is in 
terms of the existence of large gaps in the single-particle spectrum. However, 
this definition is ill-defined, the notion of a large gap being ambiguous. In 
the following, we will see that the most natural and simple interpretation 
of all the characteristic features of this figure, including the position of the 
magic numbers, is in terms of the periodic orbits of the mean-field potential. 
This simple description does not pretend to be a substitute of more elabo- 
rated many-body calculations. However, the semiclassical theory clearly does 
something very important: it grasps the essence of the physical mechanism 
responsible for the oscillatory structure of the nuclear masses, and therefore 
provides a basic tool towards a general understanding of shell effects. 



5.1 Regular motion: the spheroidal cavity 

To develop a semi-quantitative theory to explain the main features of Fig. 
a particular model needs to be specified. In the mean-field approach adopted 
here, what is required is the shape of the mean-field potential. Before select- 
ing one, we need to make some general considerations to guide our choices. In 
its ground state, the nucleus is minimizing its mass (and therefore maximizing 
23). The shell corrections are an important component of the extremization 
procedure since, locally, they control it. The type of dynamics the nucleons 
are undergoing depends on the form of the mean field potential, through the 
parameter x. The function B(A,x) has to be maximized with respect to x 
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at a fixed number of particles. Which type of dynamics is energetically more 
favorable? As we have seen in SJ2I the amplitude of the shell effects is more 
important when the motion of the nucleons is regular, by a factor yfg, as 
compared to a chaotic dynamics. To obtain a B as large, and negative, as 
possible, a mean field that produces a regular motion is therefore preferable. 

For the mean-field approach and the semiclassical techniques employed 
here to be meaningful, a large number of nucleons is required. In that case, 
experiments show that a good approximation to the mean field is a potential 
of the Wood-Saxon type, i.e., relatively flat in the interior of the nucleus 
and steep on the surface |SJ . To simplify, we will consider the potential to be 
strictly flat in the interior (and set, by convention, to zero), and infinite at 
the border. Our schematic model of the nucleus will therefore be a mean-field 
potential given by a cavity, or billiard, with hard walls: the nucleons freely 
move at the interior, and are elastically scattered off the walls. The simplest 
geometry that produces a regular motion of the nucleons, and that is of clear 
interest in nuclear physics, is a spherical cavity. We will therefore start by 
considering in some detail the mass fluctuations of a Fermi gas located at the 
interior of a 3D spherical cavity. 

Quantum mechanically, at T = the fluctuating part of the energy of a 
Fermi gas of A particles in a spherical cavity of radius R is given by, 

B sp h — — 2. €n < 1 ~ ^sph ■ (53) 

^ F i 

The normalized single-particle energies e n ,i were defined in J5J The aver- 
age part of the binding energy is given by B sp h — g s (E F /e F ){(2/15ir)e F — 
(l/8)ej + (2/9Tr)eJ ]. The sum in Eq. IJS1?|) is made over the lowest A states 
of the sphere, taking into account the spin (g s — 2) as well as the intrinsic 
{g n ,i = 21 + 1) degeneracies. The energy factor in front of the sum is written, 
for convenience, E F /e F = H 2 /2mR 2 = Eq, where 

e F = (k F R) 2 =E F /E , (54) 

E F = h k F /2m is the average Fermi energy, and m is the nucleon mass. 
Note that the normalized energy e F is equivalent, up to a constant factor, to 
the adimensional shell parameter g. Indeed, for a 3D cavity g is written, for 
A»l, 

9,Vk 2 P 

where l m in is the length of the shortest periodic orbit. For a sphere, with 
V = (4/3)7ri? 3 and £ mm = 4R, this gives, 

9 s n 0*12 9s 
9 = y(Kp-«J = -j £f ■ 

Semiclassically, the fluctuating part of the binding energy of the Fermi 
gas is given by Eqs. (|44|l and (|45J1 . The former is obtained by integrating 
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with respect to the energy up to E F (A, x) the oscillating part of the counting 
function. For a spherical cavity, AT is given by Eq. (|24(l . As usual, and to 
leading order in a semiclassical expansion, only the integration with respect 
to the rapidly varying phase factors is kept in the integral. As a result, the 
oscillatory part of the binding energy of a Fermi gas in a 3D spherical cavity 
is, 

B sph (A,x) = 4 E F g s y^ ^ -jp-sin(k F R l vw + v vw ir/2) . (55) 

m=l v =2w vw 

A vw , l vw , and u vw were defined in Eqs. l(25|) (|25|l . In A VW7 the rescaled energy 
e has to be replaced by the rescaled Fermi energy e F . 

For simplicity, we consider separately the contributions to the binding en- 
ergy of the neutrons and of the protons. The contribution of the neutrons is 
given by Eq. (|55() putting A — N. If the energy of the Fermi gas with a spher- 
ical shape is computed as a function of the Fermi energy (grand-canonical 
ensemble), it is found that the semiclassical approximation to the exact re- 
sult gives a very accurate description, of a quality comparable to what was 
obtained in JUfor the counting function (cf Fig- EJ - The aim here, however, 
is to compute B as a function of the neutron number. The dependence of 
E F (or of e F ) on N is given, according to Eq. (|45|l . by the inversion of the 
smooth part of the counting function. From Eqs. (|17J) and (J2BJ), the equation 
to invert is, 

^ 2 3/2 1 2 1/2 , . 

g s 9tt 4 37T 

When the number of nucleons is large, only the first term in the r.h.s. of 
Eq. (ISHJ may be kept. Then, 

/ q \ 1 / 3 
k F R = V^= Ij Ln ) -1-92^ 1/3 , (57) 

where we used g s = 2, the spin degeneracy. From Eq. (|54(l . the corresponding 
Fermi energy is, 

h 2 ft 2 

Ef = 7T—^ e F = 1-16 — 5- = 40 MeV , (58) 

zmR 2 totq 

where we set the nuclear radius to R = r^A 1 ^ 3 , r^ = 1.1 fori, and m — 939 
MeV. We have moreover assumed, for simplicity, an equal number of protons 
and neutrons, A = 2N = 2Z. 

Equation l|57() has a very simple and important consequence on the shell 
oscillations. The oscillating part of the binding energy is given by a sum of 
interferent terms each having, according to Eqs. (|55(l and l|57|l . a phase factor 
proportional to N 1 ^ 3 . As a function of TV, B sp h will therefore present oscil- 
lations whose wavelength grows as iV 2 / 3 . Alternatively, a plot of the nuclear 
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binding energy as a function of TV 1 / 3 , instead of N, must show constant- 
wavelength oscillations. This fact is confirmed, for the experimental data, in 
the bottom part of Fig. 

To obtain a good semiclassical description of the shell effects it is not 
enough to keep only the first term in the r.h.s. of Eq. I(56() . The full equation 
should be used to compute e F as a function of N (this is particularly impor- 
tant for the phase factors k F R in Eq. 1)55(1 . and less for the prefactors). When 
this is done, a qualitative agreement between the exact and the semiclassical 
results for the energy of a Fermi gas on a sphere is obtained (see Fig. 12 in 
Ref. 0U). The precision is however not satisfactory, and clear deviations be- 
tween the exact and semiclassical results are still observed. This is due to the 
intrinsic difficulty of computing accurately the Fermi energy as a function of 
the number of particles. The problem is particularly difficult (and probably 
one of the worst cases that can be found) for the sphere, were large degen- 
eracies are present. These degeneracies make the function E F (N,x), defined 
by inverting N — J F p(E, x)dE, ill defined. 

There are different ways to cure this problem. A first possibility is to take 
into account higher order terms in the expansion 1(42(1 in terms of //. The 
other method is to directly compute (numerically) the integral in Eq. (|40|l . 
as was done in Ref. 0J. In this way, good agreement between the quantum 
mechanical and the semiclassical results for the sphere is obtained. This is 
the method employed here to compute semiclassically the binding energies, 
shown in the figures below, as a function of N. We will continue, however, 
to refer to the less accurate (but qualitatively correct and more tractable) 
approximation l|55(l for general considerations and discussions. 

A priori, there are no adjustable parameters in the model. In particular, 
the value of k F R is determined, for a given number of neutrons, by Eq. ((56(1 . 
However, we find that a better agreement with the experimental data is 
obtained if the value of the constant in Eq. H57|l. 9tt/2 = 14.1372, is increased 
by ten percent, and changed to 15.625. The constant 97r/2 is determined by 
the geometry of the sphere. The modification can therefore be interpreted as 
a change that takes into account, (i) the effects of a soft surface instead of 
a sharp border, and (ii) the effect of spin-orbit scattering, that changes the 
lengths of the orbits in a semiclassical description 02] ■ The modification of 
this constant is equivalent to define an effective neutron number (cf Eq. (JS3), 
N e ff rs 0.91A^, and that is the way it was included in the calculations. 

We can now compare the fluctuating part of the energy obtained for a 
sphere with the experimental data. The two quantities plotted are B for 
the experimental data (as defined by Eqs.QJ and 1(51(1 '). and B sp h for the 
theory. The result is shown in Fig. [7\ It is clear that the simple model of 
a spherical cavity (with only one adjusted parameter (k F R) that has been 
slightly increased) gives a pretty good semi-quantitative description of all 
the main features of the shell effects observed in the nuclear masses: 
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Fig. 7. The fluctuating part of the mass of the 1888 nuclei of the 1995 Audi- 
Wapstra compilation (dots) (cf Fig. |SJ, compared to the fluctuating part of the 
energy of a Fermi gas on a spherical cavity (full line) . 



2) 
3) 



the period of the oscillations, that increases with N, is very well described 
by the model. The agreement between theory and experiment is particu- 
larly good for N > 25, but clear correlations are also observed at lower 
values of N. 

the spherical cavity correctly reproduces the asymmetry between sharp 
minima and rounded maxima. 

the magic numbers obtained from the spherical model, at N sp h = 2, 7, 18, 
30, 52, 83, 126, are in good agreement with the experimental values N exp = 
2, 8, 20, 28, 50, 82, 126. As was shown in [gj short orbits dominate the 
fluctuations of the energy of the gas. In order to get an analytic estimate 
of the magic numbers, we can simply consider the minima of B when only 
the triangular and square orbits are taken into account in Eq. Q55|). This 
gives the approximation, 



4) 



5) 



N„ 



(2fc+l)7T 



10.77 



9,18,33,55,84,122. 



(59) 



k is an arbitrary integer number, that has been taken equal to 3, 4, . . . , 8 
to obtain the values on the r.h.s. The constant 10.77 is computed from 
the parameters that define both orbits. 

the amplitude of the oscillations is qualitatively correct. However, in Fig.0 
only the neutron contribution is plotted. Clearly, the spherical cavity over- 
estimates the amplitude. 

The overestimate of the amplitude is larger for the maxima. This differ- 
ence is due to the well-known fact that the nucleus is deformed between 
closed shells. This is treated in more detail below. 



Figure [3 shows a direct comparison of the experimental data to the bind- 
ing energy of a Fermi gas on a spherical cavity as a function of the neutron 
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number. At each value of N, several experimental points appear in the plot, 
that correspond to different atomic elements that vary by their atomic num- 
ber Z. In its present version, the model is too simple to distinguish between 
different nuclei at a given N. The more sensible way to proceed is to simply 
compare B to the average of the experimental masses computed at a fixed N. 
This average curve is shown, with respect to the original experimental data, 
in part (a) of Fig. |H1 In part (b) the average curve obtained is compared to 
B for a spherical cavity. This figure is thus equivalent to Fig. but with 
average experimental results. 
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Fig. 8. Comparison between different theoretical approximations and experimen- 
tal values of the fluctuating part of the nuclear masses, (a) Experimental data 
(dots), with their average behavior (full line). Average experimental data (full line) 
compared to the predictions of (dashed line) (b) a spherical cavity, (c) a (possi- 
bly) deformed cavity, and (d) a (possibly) deformed cavity with a finite coherence 
length. 



The negative peaks of B define regions of more stable nuclei. In between 
magic numbers, the spherical cavity leads to neutron numbers with a positive 
contribution B. These large and positive contributions diminish the binding 
energy. A possible mechanism to avoid that effect and to increase the stability 
with respect to the spherical shape in between closed shells is deformation. 
By changing its shape, the nucleus can increase its stability by finding regions 
were B is more favorable. To complete the previous model, it is therefore nee- 
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essary to incorporate symmetry breaking degrees of freedom. In the new space 
of deformations, the mass is determined by minimizing, at a fixed particle 
number, the binding energy with respect to the deformation parameters. 

The deformation is here treated as a small perturbation of the spherical 
shape 03]. An important technical advantage of the perturbative approach 
is that it leads to a very simple semiclassical treatment, where the deforma- 
tion is incorporated as a modulation factor into the previous formulas: A p _ r is 
replaced by A p>r x M p>r (x) in Eq. (|55|) 43. x denotes the deformation param- 
eters, taken as quadrupolc (£2), octupole (X3), and hexadecapole (x^) axial 
symmetric deformations of the surface of the sphere parametrized, in spheri- 
cal polar coordinates (r, 9, </>), by r(9) = R 1 + xq + X),=2 XjPj (cos 9) . xo is 
included to ensure volume conservation, and the Pj (y) are the Legendre poly- 
nomials of order j. The explicit expressions of the modulation factors M p ^ r (x) 
for the three deformations of the periodic orbit families of a spherical cavity 
may be found in Ref. 0J . Including these factors in the semiclassical expres- 
sions, at a fixed neutron number JV, the minimum of the binding energy is 
located in the parameter space (x2,xs,Xi). 

The results of the minimization procedure for nuclei with neutron num- 
ber up to 160 are shown in part (c) of Fig. |S] With respect to a perfectly 
spherical shape, the results are greatly improved and the agreement with the 
experimental data is much better. The amplitude of the shell corrections have 
been strongly diminished. The shape and position of the negative peaks of 
B remain largely unchanged. The main difference is that the large positive 
peaks have been suppressed, and replaced by a less pronounced oscillatory 
behavior. 

As a final improvement of the model we consider inelastic processes. 
The model presented above is based on a mean field approach, were each 
nucleon (or quasiparticle) freely moves in the self-consistent mean field pro- 
duced by the other nucleons. This single-particle picture is of course an ideal- 
ized approximation of the full many body problem. In reality, quasiparticles 
have a finite lifetime associated to inelastic processes. By producing single- 
particle energy levels with a certain width, nucleon-nucleon inelastic scat- 
tering tends to wash out the coherent phenomena that lies behind the shell 
effects, and manifests as a phase-breaking mechanism. As a consequence, a 
finite phase coherence length enters now the play, that corresponds to the 
typical distance a quasiparticle travels without loosing phase coherence. 

When the finite width of the single-particle levels is included into the 
semiclassical treatment, the net output is very natural. It manifests as a 
damping factor in the sum over periodic orbits, that suppresses the contri- 
bution of long periodic orbits |3U]. This is similar to what happens to shell 
effects as the temperature of the system is raised 6 . In the approximation 
(f5K|) . the fluctuating part of the energy, including deformations as well as 
inelastic processes, is obtained by replacing A p ^ r by A PyT x M Ptr (x) x k^(£ vw ), 
where k^(£) — (£/£)/ smh(£ /£) , and £ is the phase coherent length expressed 
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in units of the radius of the sphere. For orbits whose length £ vw ^> £,, the latter 
factor produces an exponential damping of the corresponding contribution. 
In the opposite limit, when £ vw <C £, k>£ tends to one. The phase coherence 
length £ measures the typical number of times a quasiparticle can bounce 
back and forth in the mean field potential before loosing phase coherence. 

Fixing this length to £ ~ 5, the fluctuating part of the binding energy B 
is represented in part (d) of Fig. El With respect to part (c), which does not 
include decoherence phenomena, the amplitude and peaks are diminished, 
and now the theoretical description provides a very convenient description of 
the average experimental data. The agreement is less good at small neutron 
numbers N ~ 30, a foreseeable discrepancy from a theory based on a large N 
expansion. The RMS error between the theoretical curve of Fig.|H{d) and the 
average experimental data is of 1.3 MeV (taking into account only points with 
N ~ 30). This error is only a factor two larger than the best current results 
of global mass adjustments, a remarkable result considering the simplicity of 
the model. 

Before closing this section on the ability of a spheroidal cavity to repro- 
duce the experimental data, we would like to consider, within the present 
model, supershell structures. These are long-range coherent modulations 
of the amplitude of the shell oscillations. They constitute a further remark- 
able manifestation of the collective deviations of the single-particle spectrum 
with respect to its average behavior. Initially predicted for a spherical cavity 
by Balian and Bloch |30l 14*4*] , they were observed experimentally in metallic 
clusters |45] . In nuclear physics, in principle the range of variability of the 
number of nucleons is too small in order to clearly display the effect. However, 
we will see that, contrary to this expectation, there exist some indications of 
this effect in the nuclear masses. 

The simplest and more elegant formulation of supershells is based on 
semiclassical arguments. They are associated to the beating pattern produced 
by the interference of two periodic orbits of the spherical cavity that have 
similar lengths, the triangle (3,1) and the square (4,1) I^UIE]. When only 
their contributions are taken into account in Eq. I|55|l , the oscillatory part of 
the binding energy is approximated by 

B « CiV 1 ' cos I I cos I + -I , (60) 

where C is some amplitude factor that has no relevance in the present dis- 
cussion, ^ A = 3_V3 x 1.987V 1 / 3 + 9tt/4 and ip<> = 4^2 x 1.987V 1 / 3 + Htt/4 
are the phases k F R i vw + v vw w/2 associated to the triangle and the square 
orbits, respectively. This simple equation describes the supershell effect. On 
the one hand it includes the shells, given by the fast oscillations produced by 
the second cosine in Eq. I|60(l (sum of the phases). The approximate magic 
numbers Eq. I|59(l were computed by locating the minima of this term. The 
amplitude of the fast oscillations is modulated, on long range scales, by the 
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term containing the difference of the phases. These modulations produce a 
beating pattern of the shell amplitudes. The nodes of this pattern deter- 
mine the particle numbers were the shells have a minimum amplitude. From 
Eq. 1(50) 1 . the nodes are located at 



1* node 



(2n + l/2)7r 
(4\/2-3V3) 1.984 



5, 634, 3699, . 



(61) 



where n = 0, 1, 2, ... is an arbitrary integer number. Atomic nuclei of 630 
neutrons are thus necessary to observe one supershell oscillation. In spite of 
this difficulty, we can however explore if there are traces of the modulation 
in the available experimental data, having in mind that a node is predicted 
for N m 5. For that purpose, we plot in Fig. [pj the average experimental 
curve (defined in part (a) of Fig. [SJ) , and compare it to the modulation factor 
CAT 1 / 6 cos [(<p±-<p.)/2\. 
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Fig. 9. The average fluctuating part of the mass of the 1888 nuclei of the 1995 
Audi-Wapstra compilation (full line) compared to the modulation factor of the 
shell oscillations (supershell envelope, dashed line). The inset shows a closer view 
of the same curves, where the node at N sa 5 is visible. 



Although the experimental data is indeed not unambiguous due to the 
limited range of the number of neutrons, it is clear that they follow, to a 
good approximation, the supershell modulations. It should be noted that the 
increase of the amplitude observed when going from N ~ 10 to TV ~ 160 
(better displayed in the inset) cannot be explained by the overall factor N 1 / 6 
in Eq. (|60|l . It requires, as a main ingredient, the supershell modulation factor. 

The perturbative treatment of a spherical cavity, with an underlying clas- 
sical regular motion of the nucleons, allows therefore to understand in very 
simple physical terms most features of the shell fluctuations observed in the 
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atomic masses. Beyond this specific model, it is possible to use the general ar- 
guments introduced in the previous section to estimate the typical size of the 
shell effects in the nuclear masses to be expected from a generic integrable mo- 
tion. In the experimental curve, Fig. the fluctuations have approximately 
constant amplitude, with an RMS size a exp « 3 MeV. 

A theoretical prediction of the typical size of the fluctuations can be made 
from the results of $H(cf Eq. (|49|l and Ref. |7j). The variance of the binding 
energy or mass fluctuations of a Fermi gas whose corresponding classical 
dynamics is regular is given by, 

<^ es > = ^El . (62) 

We thus simply need to compute g and E c for the nucleus. The length of the 
shortest periodic orbit can be estimated as twice the diameter of the nucleus. 
Its period is T m j n = AR/v F , where R is the nuclear radius and v F the Fermi 
velocity. Then, according to the definition (f 1 1|) . 

h hv F ttE f 



T min 4R k F R 

From Eq. IJ57J) (or its generalization to an arbitrary shape), with N = A/2, 
we have k F R « l.hA 1 / 3 . Therefore, using moreover Eq. JSSJ, we obtain, 

80 
E * = 37S MeV • ( 63 ) 

On the other hand, from Eqs. (|15|l and JTSJ, putting A/ = A, it follows that 



and hence, 

9 



^ « 3 x A 2 '* . (64) 

o 



In Eq. (|62|l . the A dependence between E^ and g cancels exactly, and the 
RMS of the mass fluctuations of a regular motion is, finally, 



(B? eg ) « 2.9 MeV , (65) 

which is in good agreement with the experimental result. For a generic reg- 
ular dynamics, the size of the mass fluctuations is therefore expected to be 
constant, e.g. with no dependence on the number of particles. In this respect, 
the spherical cavity is special. In fact, Eq. 125J I shows that the amplitude of 
all the periodic orbits (except the v — 2w family, which is of lower weight) 

scales as e F ~ A 1 / 6 (cf Eq. (|57|l ). Therefore, the size of the fluctuations for a 
spherical cavity scales as a S ph ~ A 1 ' 6 . In atomic nuclei, this prediction is not 
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generically valid, due to deformations. But this scaling is expected to hold for 
the strength of the peaks associated to the magic numbers, where the nucleus 
is well described by a spherical shape. From a dynamical point of view, the 
specificity of the sphere compared to a generic regular motion comes from 
the fact that the codimension of the v > 2w families of periodic orbits of the 
sphere is three, whereas the generic codimension for 3D integrable systems is 
two. 

5.2 Chaotic motion 

The description of the mass fluctuations in terms of a regular nucleonic mo- 
tion is globally quite satisfactory. However, the semiclassical methods em- 
ployed here only provide a semi-quantitative description. Actual calculations 
of the masses are in fact much more sophisticated. Different methods, more 
or less phenomenological, have been used to produce global systematics of 
the nuclear masses. These are based on shell model approaches, mean field 
theories (with or without correlations, using different effective forces), or 
phenomenological liquid drop models plus shell corrections. The number of 
adjustable parameters in these theories is considerable, of the order of 20-30. 
See Ref. 02] f° r a recent review that describes the different approaches. In 
spite of the sometimes very different conceptual basis of these models, there 
is a peculiar common feature that emerges. Over the whole mass chart of 
known nuclides, the different models possess a comparable precision, of the 
order of 400 — 700 keV. A typical RMS error with respect to the experimental 
data, taken from Ref. 07J> is shown in Fig. E3 (dots). An error of the order 
of 500 keV, which decreases with A, is observed. 

Is there a way to understand this deviation? The fact that different models 
have approximately the same precision points to the fact that there may be 
a fundamental physical mechanism at work here. In the present framework 
based on a semiclassical description of the motion of the nucleons, a natural 
scheme to explore is the presence of some chaotic motion. Although we have 
argued that a regular motion of the nucleons is energetically more favorable 
than a chaotic one, there are several reasons that make the presence of the 
latter inevitable. The first one is that a perfectly regular motion is rather 
unusual in nature (it holds for the two-body problem with central forces, but 
the three body is already not integrable). A mixed dynamics, where regular 
and chaotic motion coexist, is generically expected to occur. There is no 
simple reason to believe that the full many-body nuclear problem, with of 
the order of 100 interacting particles, will not follow the generic rule. The 
second argument in favor of the presence of chaos is also fundamental. As 
we will argue below, there are indications that what is really computed by 
assuming some chaoticity in the nucleon's motion are correlation effects that 
are beyond the mean-field description of the nucleus. 

We therefore assume that the phase-space of the nucleons at Fermi en- 
ergy is dominated by regular components but contains, nevertheless, some 
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Fig. 10. The RMS of the difference between computed and observed masses as a 
function of the mass number (dots, from Ref. |47)1. compared to the RMS of the 
mass fluctuations due to a chaotic motion. Taken from Ref. 0. 



chaotic layers. The idea is to compute the contribution to the mass of these 
chaotic layers, and then compare it to the deviation observed between the 
experimental data and the best theoretical predictions. From a semiclassical 
point of view, the most simple approximation that can be made in the generic 
case of a mixed dynamics is to split the sum over periodic orbits in Eq. I|44J) 
into two terms, one from the regular part, the other from the chaotic phase 
space sector. The mass fluctuations are now written as, 

B = B reg + B ch . (66) 

The two terms on the r.h.s. are, from a statistical point of view, independent, 
(B re gB c h) = 0. This happens because, as already pointed out in jJH the 
dominant contribution to the energy comes from the short orbits, with t p <C 
t h . Since the orbits contributing to each term are different because they lie 
in different phase-space regions, the cross product of two different orbits 
vanishes by the averaging procedure (assuming that the actions of the orbits 
are incommensurable). 

We now need to evaluate the variance of the shell corrections that orig- 
inate in the chaotic layers of a nucleus, given by the term B c h in Eq. (|66|l . 
We want to make a statistical estimate which, as discussed before, has the 
advantage of requiring a minimum amount of information about the system. 
This can be done from the general result obtain in 2] Eq. Ij49(l with (3 = 1 
(assuming time reversal symmetry), 

(%h) = ^K ■ (67) 
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An estimate of E c is now required. In fact, the result obtained previously, 
Eq. (|63fl . can be used here. The reason is that E c is defined as the energy 
conjugate to the shortest periodic orbit. In a system were regular motion 
coexists with chaos, the shortest orbit lying in the regular or in the chaotic 
components will have approximately the same period (or length) at fixed en- 
ergy E F (A, x). Therefore, the scale E c associated to each of the contributions 
in Eq. Hfifijl will be similar. 

No other information is needed to estimate the size of the contributions 
to the mass that originate from the chaotic layers. Using Eq. (|63fl . we obtain, 



<rch = \/(B 2 ch )*-^MeV. (68) 

This result has to be compared with the RMS of the difference between the 
experimental and the computed masses, which we take, as mentioned before, 
from Ref. 07J- The comparison is shown in Fig. 1101 

The agreement is quite good. The amplitude in Eq. (|68(l is uncertain up to 
an overall factor of say, 2. It can be varied by increasing slightly the period of 
the shortest orbit (we have chosen twice the diameter of the nucleus, which 
is the shortest possible length; any modification, assuming for example a 
triangular-like shape, or any other shape, will increase the length and period, 
and therefore will diminish a c h), and by a more appropriate inclusion of spin 
and isospin (this increases a c h by a factor 2 if these components are treated as 
uncorrelated). The A dependence of the error is very well fitted in the region 
A ~ 75, with deviations observed at lower mass numbers. This is consistent 
with the expected loss of accuracy of the semiclassical theories for light nuclei. 

There are several features of the present theory - that interprets the mass 
deviations in terms of chaotic motion - that make its predictions reliable. 
First of all, Eq. (|6Tjl contains only one physical parameter, the period of 
the shortest chaotic periodic orbit, which is a function of A because the size 
of the nucleus increases with the mass number. It has no dependence on 
the relative size of the chaotic region, i.e., on the fraction of phase space 
occupied by chaotic motion. Without this quite remarkable and important 
fact we would have been forced to estimate that fraction, something that is 
hardly possible with present knowledge, despite the efforts in this direction 
since the pioneering work of Ref. [48|. Mathematically, this is due to the 
fact that, in chaotic systems, the dominant part of the form factor in the 
integral (|4"TI) is independent of the Heisenberg time (cf Eq. JSJl). But it is 
precisely that time that has information on the corresponding level density 
6 = ~JT X (the total average level density of a mixed system is the sum of 
the average level densities of the regular and chaotic components. Each of 
these densities is proportional to the phase space volume occupied by the 
corresponding motion). Therefore, for the chaotic component the integral 
(|47(l has no information on the relative phase space volume occupied by 
the chaotic layers. This information enters only as a correction, through the 
parameter g in Eq. (j49(l . 
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Equation I|tj7|) is in fact quite robust. It is not only independent of the 
chaotic phase-space volume, but it is also valid for arbitrary dimensions. Af- 
ter a system-dependent transient (schematically represented in Fig. [3J, the 
linear growth of the form factor in Eq. 147|) for chaotic systems with /3 = 1 is a 
quite general feature, independent of any information of the system. This fact 
suggests a very natural origin of these fluctuations. Although our analysis is 
based on a single-particle picture, it can be extended to the full many-body 
phase space (of dimension 6xA). In that space, and in a semiclassical picture, 
to first approximation the point representing the system follows a very simple 
trajectory, driven by the regular mean-field that dominates the motion. On 
top of that, it is likely that the residual interactions, not taken into account 
in that approximation, induce chaotic motion. The presence of chaotic orbits 
would then introduce additional long-range modulations. Eq. Ij67(l would then 
still be valid to evaluate the amplitude of those modulations, with T min the 
period of the shortest chaotic orbit in the multidimensional space. Rough 
estimates indicate that T m i n is comparable to the three-dimensional non- 
interacting period (this corresponds, for example, to collective and coherent 
oscillations of the particles) . Therefore Eq. (jf>T|> presumably gives an estimate 
of the mass fluctuations arising from neglected many-body effects. This de- 
serves further investigation, though some work in this direction has already 
been undertaken 49 . 



6 Conclusions 

We have presented in these lectures a study of shell effects on nuclear ground 
states within a unified framework, namely periodic orbit theory. In a mean- 
field approximation of the nuclear dynamics, the single-particle motion is 
quantized. The regular or chaotic nature of the classical single-particle mo- 
tion imprints the single-particle spectrum in two different ways. The first one 
occurs at the scale of the mean level spacing S, and produces different (but 
universal) statistical fluctuations. The second acts on a much larger scale 
E c , set by the inverse of the time of flight across the system. A bunching of 
the single-particle levels produced by the short periodic orbits is observed 
on this scale, whose intensity depends on the regular or chaotic nature of 
the motion. These two features of the single-particle spectrum influence the 
probability distribution of the quantum fluctuations of thermodynamic func- 
tions of the Fermi gas in different ways. In the particular case of the mass, 
the fluctuations are controlled by the short periodic orbits. 

Without need of heavy and numerically expensive many-body calcula- 
tions, the (phenomenological) semiclassical theory, although less precise, pro- 
vides a good semi-quantitative explanation of most of the features of the 
experimental data on nuclear masses. A simple model of regular classical 
motion based on a spherical cavity (including the possibility of deformations 
and decoherence effects) gives already a good agreement. Within this model, 
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the average behavior of the experimental data is reproduced, for nuclei with 
more than 30 neutrons, with a RMS error of 1.3 MeV. This error is only a 
factor two larger than the best current results of global mass adjustments, a 
remarkable result considering the simplicity of the model. This stresses that 
the main physical effect behind shell effects is purely kinematic, with spin- 
orbit or other interactions having a small overall influence. The latter were 
phenomenologically taken into account by adjusting the value of k F R. Within 
this simplified model, long range modulations of the shell amplitudes (super- 
shells) were also investigated. Although, due to a limited particle number, this 
effects cannot unambiguously be displayed, clearly the present experimental 
data show some indications of supershell structure in the atomic masses. 

When the experimental data is compared with the best global theoretical 
calculations over the whole mass spectrum, deviations of the order of 500 keV, 
which decrease with A, are observed. We interpret this deviation as due to 
the presence of some chaotic components in the phase space dynamics of the 
nucleons. Our computations require very little information on the system. In 
spite of that, they are in good agreement with the deviations observed. Our 
interpretation does not imply the existence of a "chaotic" insurmountable 
obstacle that, intrinsically, limitates our understanding and predictive power 
to compute the nuclear masses (see, in this respect, Ref. |SS1)- It simply says 
that at these scales (of the order of 500 keV) a new effect appears that needs 
to be taken into account properly. The picture we are suggesting is coexistence 
of order and chaos as produced, for instance, by residual interactions. Further 
evidence should be given to confirm this scenario. In particular, the theory 
predicts autocorrelations in the total energies 6^ , as well as the effect of the 
presence of chaotic layers on the nuclear level density (as a function of the 
excitation energy). 
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